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Abstract. The UH8K wide field camera of the CFHT 
was used to image 0.68 deg 2 of sky. From these images, 
^20,000 galaxies were detected to completeness magni- 
tudes V < 23.5 and / < 22.5. The angular correla- 



galaxies as a function of their redshifts, spectral types and 
star formation rates (Kauffmann et al., 1999a; 1999b). By 
measuring redshifts for 10 5 or more galaxies, the next gen- 
eration redshift surveys such as the VI RMOS survey (pc 



tion function of these galaxies is well represented by the Fevre et al., 1998), the DEEP survey (Davis and Faber 



parameterization uj{9) — Aw 9 ■ The slope 5 ~ —0.8 1998), and the LZT survey (Hickson et al., 1998; see also 



shows no significant variation over the range of magni- 
tude. The amplitude Aw decreases with increasing mag- 
nitude in a way that is most compatible with a ACDM 
model (Oq = 0.2, A = 0.8) with a hierarchical clustering 
evolution parameter e > 0. We infer a best-fit spatial cor- 
relation length of r o — 5.85 ± 0.5 ft. -1 Mpc at z = 0. 
The peak redshift of the survey (/ < 22.5) is estimated 
to be Zp ea fc ~ 0.58, using the blue-evolving luminosity 
function from the CFRS and the flat A cosmology, and 
r o( z peak) — 3.5 ± 0.5 hr 1 Mpc. We also detect a signifi- 
cant difference in clustering amplitude for the red and blue 
galaxies, quantitatively measured by correlation lengths of 
r 0Q = 5.3 ± 0.5 h' 1 Mpc and r 00 = 1.9 ± 0.9 h" 1 Mpc re- 
spectively, at z — 0. 

Key words: galaxies - angular correlation function - clus- 
tering - survey 



1. Introduction 

For the past decade, the study of the spatial Large-Scale 
Structure (LSS) of the universe has become an major tool 
for constraining the cosmological models. In particular, 
provided many assumptions on how morphological type 
correlates with colour, how mass is correlated with opti- 
cal luminosity, and how local density correlates with mor- 
phology, recent CDM hierarchical N-body simulations and 
semi-analytic models of galaxy formation are able to make 
tentative predictions on the clustering evolution of the 
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http://www.astro.ubc.ca/LMT/lzt/index.htm]) will allow 
detailed studies of the large-scale clustering and its evolu- 
tion tO Z r-J 1. 

Until these surveys are completed, the measurement 
of the 2-point angular correlation functions to (9) of large 
photometric galaxy samples remains the best alternative 
to constrain galaxy clustering at z > 1. The major caveat 
of u>(0), as opposed to the 2-point spatial correlation func- 
tion £(r), is that it probes the projection of a 3D distribu- 
tion of the galaxies onto the 2D celestial sphere, i.e. one 
can not tell whether a given galaxy is a faint nearby object 
or a bright remote one. As a consequence, u>(9) is sensitive 
to the effects of both the intrinsic 3D clustering and the 
luminosity evolution (LE) of galaxies for a given set of cos- 
mological parameters. To avoid this degeneracy, one must 
choose between two approaches to extract sensible infor- 
mation from u>(6). First, one may use past observations 
to assume a scenario of galaxy evolution with given LEs 
and redshift distributions for each galaxy population, and 
then deduce the clustering evolution. A second approach 
would be to assume a clustering scenario, cosmological pa- 
rameters, and to measure the evolution of the correlation 
function in order to validate the theoretical LE used to 
model the galaxy counts, e.g. Roche et al. (1993). In this 



paper, we favour the first approach. We use the Canada- 
France Redshift Survey (CFRS; Lilly et al. 1996) luminos- 
ity function and redshift distribution to invert the angular 
correlation function with Limber's formula (cf section^ on 
modeling of lu[9], £[r] and r ) to compute the spatial cor- 
relation length from the amplitude of lo(9). This approach 
has several limits which are discussed in section |7]. 

An extensive literature covers the evolution of cluster- 
ing using to (9). The first attempts to d escrib e clustering 
using counts in cells is due to Limber ( 1954 ). The two- 
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point correlation function as a statistical diagnosis of clus- 
tering has been popularized in astrophysics by Hauser & 
Peebles ( 1973|) and app lied to the Zwicky Catalog (Pee 



bles ai ^d Hauser, 1974 ). Since these pioneering studies 



the method, fully described by Peebles (198C), has been 
applied to many photographic catalogues in diverse pho- 



tometric bands e.g. see Groth & Peebles (1977), Koo & 
Szalay (|1984[), Maddox et al. fll990b|) where the clustering 



of local galax ies is studied on large angular scales. Us- 
ing Limber's ( 1953 ) formula relating tu(9) to the spatial 
correlation function £(r), these studies establish that the 
spatial clustering of local galaxies can be parameterized 
as a power law, £(r) = (r/ro)™ 7 , where the correlation 
length r ~ 4 - 8/i _1 Mpc (h = H /100) and the slope 
7-1.7- 1.8. 

A second generation of studies based on small-scale 



et al, 
terson 



CCDs, probes smaller areas to deeper magnitudes (Ef- 
stathiou et al., 1991 ; |Campos et al., 1995 ; Ncuschaefer 



1992; Neuschaefer et al., 1995; Roukema and Pc- 
1994; Braincrd et al., 1995|; Braincrd and Smail. 



1998|~p jraincrd et al., 1999[ [Hudon and Lilly, 199C|; [Lid 



1996; 'Voods and Fahlman, 1997), therefore allowing to 



man and Peterson, 199(3); [Roche ct al., 1993| ; [Roche et al 



measure the evolution of the correlation length to red- 
shifts z < 1). The most recent studies take advantage of 



large area CCD detectors ( Roche and Eales, 1999 ; Post 
man et al., 1998] ) to measure the angular corre lation func 



tion, and of the use of photometric redshifts ( Koo, 1999 ) 
to estimate the spatial correlation function from deep pho- 



tometric surveys ( Vlllumscn ct al., 1997 ; Arnouts et al 
1999fj~|Despite these numerous studies, our knowledge of 
the clustering of galaxies is still rudimentary. The main 
trends are that while a mild luminosity evolution seems 
to be required to explain faint number counts in the I- 
band, weak or no evolution in the galaxy clustering with 
redshift is detected out to z ~ 1. 

The existing surveys measuring the galaxy angular cor- 
relation function, and covering a broad range of magnitude 
bands and limits, constrain the value of r {z pea k) to the 
range 1.5 — 4.5/i~ 1 Mpc for 0.5 < z peak < 1. The disper- 
sion is mainly due to the uncertainty in our knowledge 
of the luminosity functions and redshifts distributions for 
the different galaxy types at z <~ 1 (section 0) , and possi- 
bly to the varying selection biases from survey to survey. 
To be consistent with the values of roo measured from 
the nearby redshift surveys and ranging from 4/i~ 1 Mpc 



to 8/r 


^Mpc (dc Lap 


□arent et al., 1988; Loveday et al., 


1992 
1998 


; 


Jole et al., 1994 




Tucker et al., 1997; 


Ratcliffe et al., 


; 


luzzo et al., 1998 


) , most observations of the galaxy 



clustering favour either constant or increasing clustering 
with time in proper coordinates, which is consistent with 



N-body simulations of CDM hierarchical universes ( Davis 
et al., |L985| ; |Baugh ct al., 1999| ; [Hudon and Lilly, 1996[ ) 



Note that only a few redshift surveys allow a direct 
study of the evolution of the spatial clustering (Lilly 
et al., |1995| ; |Connolly et al., 1998[ |Arnouts et al", 1999| ; 



jCarlberg et al., 200C ): these survey measure ro(z pea k) ~ 
1.4 - 4.5 ft^Mpc for 0.5 < z peak < 1. We u nderline that 
except for the results of Carlberg et al. ( [2000 ), correspond- 
ing to the high value of ro(z pea k), the limited area of the 
mentioned surveys make them very sensitive to cosmic 
variance, and the corresponding results on the correlation 
function must be taken with caution. 

Moreover, the existing analyses have not yet answered 
convincingly to the following two questions: Is there an 
evolution of the angular correlation function slope S = 
7 — 1 at faint limiting magnitudes? And do red-selected 
objects and blue-selected objects show a true difference in 
3-D clustering? In addition to providing another measure 
of the galaxy clustering at z ~ 0.5, the new sample pre- 
sented here allows us to address these questions. The pa- 
per is organized as follows, the observations are described 
in section ||, section || presents the data reduction, sec- 
tion H addresses the star/galaxy separations and counts. 
Section Jq details the analysis of the correlation function, 
section ^[gives the results, and section ^ provides a discus- 
sion of our results and a comparison with previous work. 

2. Observations 

The data were obtained during the spring of 1998 using 
the prime- focus wide-field UH8K CCD Camera (Metzgei 
|ct al., 1995 ) at the Canada-France-Hawaii 3.6 m Tele- 



scope. The camera is a mosaic of 8 2kx4k CCD chips cov- 
ering a total area of ~ 28' x 28' with a scale of 0.206"/pixel. 
We observed four fields with R.A. offset of 23', 23', and 17', 
hence overlapping by 5', 5' and 10' (Table [l]). The obser- 
vations were done in bands V (Johnson) and I (Cousins). 
The total area of the survey is 0.68 deg 2 . The limiting 
magnitudes are V < 24 and / < 23 (see section |3.3| on 
the completeness). The central galactic coordinates are 
/ = 126.2°, b = 68.2°. At such high g alactic latitudes, the 
redde ning is E(B — V) < 0.03 mag ( Burstein and Heilcs 



1982 ). This is an upper limit because we observed in V 
and I (reddenings A -1 ). We neglected both the absolute 
reddening and the relative reddening between the fields. 
Table gives the characteristics of each field. The seeing 
is between 0.7" and 0.8" for both filters during the whole 
night. The exposure time is 1200 sec for all fields. Fig- 
ure [| shows the map of the ~19,500 galaxies detected to 
/ < 22.5 (limiting magnitude of the correlation analysis, 
cf section |3.3| ). In spite of the bad cosmetics of CCD#4 
(visible in the upper right corner of Fig. |l|) , the 5' overlaps 
provide homogeneous sampling of the area, except for the 
regions containing bright stars (empty circle in the up- 
per middle) and the gaps within the UH8K Mosaic CCD 
chips. 

3. Data reduction 

The data were reduced using FLIPS by J.C. Cuillandre as 
part of the calibration service of the 1998 season of UH8K 
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Fig. 1. Map of 19,506 galaxies to I < 22.5. Bad cosmetics, intervals between CCD chips, and saturation artifacts due 
to bright stars are masked (shown as empty areas). 



Table 1. Field characteristics 



Field 


02000 


#2000 


airmass 


airmass 




hr min sec 


o / // 


in V 


in I 


UD02 


12 37 55.3 


+49 08 50.6 


1.491 


1.161 


UD03 


12 40 15.6 


+49 08 52.3 


1.371 


1.148 


UD04 


12 42 35.8 


+49 08 54.1 


1.285 


1.154 


UD05 


12 44 06.3 


+49 08 55.3 


1.224 


1.176 



Camera, Cuillandre fll998a[ |l998b| ). FLIPS suppresses the 
dark and bias using the CCD overs-cans and flattens the 
response of the 8 detectors using a flat-field made by com- 
bining ~ 300 images of 18 nights of observation with the 
UH8K Camera. The final images have a corrected sky flux 
showing variations of less than 1% within each CCD im- 
age. The photometric standards were pre-reduced follow- 
ing exactly the same steps. The final dataset consists of 64 
2kx4k frames (32 in V, 32 in J), and a set of photometric 
standard stars. 

P hotometry was performed using the SExtractor Pack- 
age ( Bertin and Arnouts, 1996 ) which provides Kron-like 
elliptical aperture and isophotal fluxes, (X, Y) coordi- 
nates, position, elongation and stellarity class for all ob- 
jects above a given threshold (we choose a threshold of 
25 contiguous pixels above 1.5 a of sky value in V and 
/) . The resulting SExtractor files were then calibrated for 
astrometry and photometry. 

3.1. astrometry 



coordinates of most objects in our fields to a red mag< 20 
with accuracies of < 0.5", and because it is easily ac- 
cessible via on-line astronomical databases. The astrom- 
etry is done separately on the 64 frames, with IRAF ge- 
omap/geotran second-order Legendre polynomials. A ra- 
dial correction is applied previously on the {X, Y) coordi- 
nates to correct the prime-focus optical corrector distor- 
tion, whose equations were provided by J.C. Cuillandre 
( 1996p . If R is the actual radial coordinate of an object 
from the center of the UH8K mosaic (in mm) and r its 
observed radial coordinate (in mm), then the shift r — R 
due to the corrector is 



r-R = 9.07 x 10" 7 r 3 + 2.06 x 10" 12 r 5 
r = 0.07284 0(1 + 2.593 x 10~ 9 6> 2 093 ). 



(1) 



The USNO-A2.0 Astrometric Catalogue flMonet, 1998|) is 
used as the astrometric reference as it gives the equatorial 



9 is the angular distance of the object from the center in 
arcsec. According to cq. (0), the radial distortion is 0.13 
pixel at a radius of 3', 0.32 pixel at 4', and it becomes 
non-negligible for radii greater than 6' where the distor- 
tion is greater then one pixel (e.g. the radial distortion at 
14' is 14.1 pixels). The overlaps between the fields allow 
us to verify the importance of the optical distortion cor- 
rection on the accuracy of the final astrometry. Figures 
and H show the errors Aa and Ad in X and Y directions 
for the overlap between the fields UD03 (CCD #0) and 
UD04 (CCD #6) without correction of distortion (Fig. [§) 
and with correction (Fig. ||). One can see that both the 
systematic errors and the random errors are bellow 0.5" 
after correction. We are not able to remove completely the 
systematic shifts between the fields (Fig. || frame a and b). 
This may be due to unknown misalignments between the 
chips of the mosaic. One can only minimize the system- 
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Fig. 2. Differences between the equatorial coordinates of 
the objects of the overlapping region of the fields UD03- 
CCD#0 and UD04-CCD#6 (in arcsec) if no correction is 
applied for the prime focus optical corrector distortion, 
(a) shows the right ascension difference Aa vs the CCD X 
axis, (b) shows Aa vs the CCD Y axis, (c) shows the dec- 
lination difference AS vs the CCD X axis, and (d) shows 
AS vs the CCD Y axis. Important systematic effects are 
seen in frame (c) and (d). 



atic effects to sub-arcsecond values. The final equatorial 
coordinates are precessed to J2000. 

McCracken (priv. comm. 2000) used a finer and more 
complex approach with the Canada-France Deep Field 
(CFDF) (~ 1 deg 2 covering the CFRS in UBVRI, the 
fields in BVRI bands were observed with the UH8K Cam- 
era) by establishing an internal Word-Coordinate System 
to which all pointing are referred. The claimed dispersion 
is 0.3 pix = 0.06" over the entire camera. 

3.2. photometry 

A rigorous photometric calibration is usually performed 
by applying a linear transformation on the instrumental 



Fig. 3. Same as Fig. g with the correction of the prime 
focus optical corrector distortion of eq. ([!]). 

(observed) magnitude to of the objects (given by Sextrac- 
tor Kron-like elliptical aperture fluxes), 



M 



A m x X + ^rn x Am + TOq, 



(2) 



where M is the standard magnitude, A rn is the extinction 
coefficient, \ 1S the airmass, k m is the true-colour coeffi- 
cient, Am is the true colour and too is the zero-point. 

In principle, the coefficient A m must be derived from 
the observation of a field of standard stars at three differ- 
ent airmasses, and the coefficients k m , from a large range 
of star colours. Because the zero-points too and the colour 
coefficients may vary from CCD chip to CCD chip, one 
should also observe the standard field separately in the 8 
CCDs of the mosaic at least three times through the night. 

For the UH8K Camera, this task is clearly beyond the 
observer's reach because the reading time of the UH8K 
camera is too long and the whole night would not be suffi- 
cient to observe the fields required for a proper photomet- 
ric calibration. In fact, the dataset provided by the service 
observing consists of one field of standard stars SA104 



(Landolt, 1992), observed in the middle of the night, in / 



and V. This minimal observation only allows one to mea- 
sure an average zero-point over the whole mosaic in each 
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filter by combining all the standard stars (one or two per 
CCD chips). 

We derive the colour coefficients ky and ki and the 
magnitude zero-points Vq and Iq by least-square fit of the 
colour transformation equations to the SA104 sequence: 



(V- 
(I 



my) 
- mi) 



kyx(V- I) 

ki x (V -I)- 



-V 



(3) 
(4) 



The airmass correction is included in the my , to/ mag- 
nitudes, and we use the standard CFHT values for a 
thick CCD, Ay = -0.12 and Aj = -0.05 (http: 
//www. cfht.hawaii.edu /Instruments /Imaging /FOCAM 
/appen.html^F) The results appear in Table [| For com- 
parison, the Table also lists another UH8K measurement 
obtained by J.C. Cuillandre as part of the calibration ser- 
vice of the UH8K 1996 season; no corresponding errors are 
provided. 

For further control, we also derive the ky and ki coef- 
ficients from a synthetic standard sequence: synthetic "in- 
strumental" magnitudes are computed for a set of stellar 
template spectra of different colours by a routine which 
uses the theoretical transmission curve of the instrument 
optics + filter (Arnouts, priv. comm. 1999). Application 
of eq. (3) and (4) to the synthetic sequence yields the 
"synthetic" ky and ki, also listed in Table [|. The corre- 
sponding zero-points Vq and Iq are derived in a second 
step by application of eq. (3) and (4) to all stars of the 
SA104 sequence, this time with the "synthetic" ky and 
kj. 

The colour coefficient ky derived from the synthetic 
sequence is compatible with the other values listed in 
Table |^. A significant dispersion appears in the fc/ mea- 
surements, probably because this coefficient is small. The 
zero-points Vq and Iq listed in Table || are also consis- 
tent within the error bars and are poorly dependent on 
the colour coefficients. Because the colour coefficients ky 
and kj derived from SA104 display large errors, we choose 
to adopt the colour coefficients from the synthetic se- 
quence and the corresponding zero-points calibrated on 
SA104: ky = -0.035 ± 0.001 and ki = 0.005 ± 0.0003; 
V = 24.58 ± 0.018 and 7 = 24.76 ± 0.014. 

For all observed objects in the catalogue, the instru- 
mental magnitudes are converted into standard V and I 
magnitudes by re- writing the colour equations eq. (3) and 
(4), written in terms of the standard colour V — I, into 
functions of the measured instrumental colours (my —mi): 



, ky x (my ~ mi) 
V = m v -Ayx x + -, \ r " + v o (5) 

1 + ki — Ky 

T A , kl X ( mv _ mi ) , T fa\ 

I = mi - Ai x x + — ttt T, — " + 7 o ■ 6 

1 + Ki — ky 

mi , my are the observed magnitudes (Kron elliptical aper- 
tures), and ky fe/ V and J are the values labeled "syn- 
thetic" in Table 0. The other coefficients have the same 



meaning as those of eq. (Q), and we use, as in eq. (3) 
and (4), the standard CFH values Ay = —0.12 and 
Ai = -0.05. 

The astronomer in charge of the service observing 
stated that the night was clear with only thin cirrus visi- 
ble near the horizon at sunrise (Picat, priv. comm. 1998). 
The overlapping regions between the various mosaic fields 
observed can be used to estimate the possible variations 
in the zero-points during the night. The fields were ob- 
served in the following sequence: UD02-I, UD03-I, UD04- 
I, UD05-I, UD05-V, UD04-V, UD03-V, UD02-V. Figure | 
plots the average of the magnitude differences Amag for 
the bright objects (with / < 19 and V < 20) detected in 
the overlapping CCD regions as a function of sidereal time. 
For each mosaic field observed, there are 2 to 4 CCD's pre- 
senting an overlap with a CCD within another field, and 
each average Amag is plotted at the sidereal time of the 
first observed overlap. 

A small systematic variation of the average Amag, de- 
noted < Amag>, with sidereal time is detected in the I 
filter, and possibly in the V filter. A field-to-field correc- 
tion of the zero-points is done to account for these small 
variations during the night: we apply to each field a cor- 
rection in its zero-point measured by the < Amag> in Fig. 
[| taking the fields UD02-I and UD02-V as references, and 
following the sidereal sequence; the same zero-point cor- 
rection is applied for all the CCD of each mosaic field. 
The residual magnitude variations measured after correc- 
tion in the CCD overlaps are a ~ 0.05 mag in both the V 
and I filters. These put an upper limit on the variations in 
the zero-points and colour coefficients between the differ- 
ent CCD's of the mosaic which are not accounted for in 
the present analysis. Note that a gradient remains in the 
V data (see section 6.3), which will be removed using the 
variation of the average galaxy number counts variations 
with right ascension. 

To evaluate the final photometric errors in the ob- 
tained catalogue, the V and I magnitudes of objects in 
the overlapping sections are also compared individually. 
Figure || gives the residuals in the V and I bands versus 
magnitude for all objects having V-I colours (cf section 
3.3) and Table || gives the corresponding standard devia- 
tions. The V — I errors are taken to be ~ \/2 ctav- The 
large dispersion at bright magnitudes in Fig. || is due to 
saturation effects. The "tilted" variation of the residuals 
with magnitude for residuals larger than 0.5 mag in abso- 
lute values (in both the V and I bands) is caused by the 
following effect: each object in the overlaps is given the 
magnitude measured in one of the overlap, arbitrarily. If 
the average of the 2 magnitudes in the overlap were used, 
this "tilt" would vanish. However, this affects only ~ 0.5% 
of the objects, and we consider that it would make no sig- 
nificant difference in any of the results reported here. 
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Table 2. Colour coefficients k for the V and / bands 
derived from the standard field SA104, J.C. Cuillandre's 
calibration service, and S. Arnouts' routine (labeled Syn- 
thetic, cf text). The average zero-points <zero-point> 
were derived using the associated k. 



Source 


Filter 


k 


<zero-point> 


SA104 


V 


0.0 ±0.05 


24.62 ± 0.06 


SA104 


I 


0.041 ± 0.041 


24.72 ± 0.05 


Cuillandre 


V 


-0.033 


24.629 


Cuillandre 


I 


-0.066 


24.721 


Synthetic 


V 


-0.035 ± 0.001 


24.58 ±0.018 


Synthetic 


I 


0.005 ± 0.0003 


24.76 ± 0.014 



d 



E o 

IN 

d 



12 



1 4 



16 



Sidereal Time (hours) 
Fig. 4. Average magnitude differences Amag of the bright 
objects present in the overlapping regions versus the side- 
real time. The /-band data are shown as filled symbols, 
and the y-band data as open symbols. Time variations 
are always smaller than chip-to-chip zero-point errors. 





Fig. 5. Photometric residuals AV^ and AI in the overlap- 
ping frames versus magnitude V and /. The vertical lines 
indicate the completeness limits. The associated standard 
deviations are given in Table |. 

3.3. Magnitude and colour completeness 

Three catalogues are generated from the data: one cat- 
alogue in the V band, one in the I band, and one cata- 
logue containing objects with measured V — I obtained by 
merging the two catalogues. All objects whose centroids 



Table 3. Photometric errors in V and / (see Fig. ||). 
The errors combine random noise, residual variations from 
field-to-field in the zero-points during the night, and the 
uncorrected variations of the zero-points and colour coef- 
ficients from CCD-to-CCD within the mosaic. 



V 


N 




<AV> 




N 


i^a; 


< AI > 


15.5 


2 


0.018 


0.029 


15.0 


3 


0.030 


-0.017 


16.0 


14 


0.061 


-0.021 


15.5 


23 


0.046 


0.037 


16.5 


14 


0.076 


-0.036 


16.0 


23 


0.061 


0.006 


17.0 


16 


0.067 


-0.002 


16.5 


37 


0.016 


0.023 


17.5 


27 


0.048 


-0.013 


17.0 


45 


0.030 


0.014 


18.0 


33 


0.045 


0.019 


17.5 


51 


0.071 


0.025 


18.5 


39 


0.042 


0.010 


18.0 


103 


0.027 


-0.005 


19.0 


69 


0.048 


0.002 


18.5 


108 


0.047 


0.011 


19.5 


106 


0.060 


-0.010 


19.0 


203 


0.097 


-0.029 


20.0 


138 


0.105 


-0.000 


19.5 


340 


0.052 


0.015 


20.5 


201 


0.120 


-0.011 


20.0 


498 


0.122 


0.043 


21.0 


321 


0.143 


-0.017 


20.5 


752 


0.101 


0.036 


21.5 


502 


0.188 


-0.030 


21.0 


1077 


0.138 


0.063 


22.0 


856 


0.233 


-0.039 


21.5 


1360 


0.157 


0.075 


22.5 


1244 


0.274 


-0.032 


22.0 


1666 


0.238 


0.072 


23.0 


2005 


0.333 


-0.049 


22.5 


1921 


0.307 


0.091 


23.5 


2596 


0.379 


-0.015 


23.0 


1712 


0.407 


0.120 


24.0 


2247 


0.471 


0.183 


23.5 


1003 


0.505 


0.138 


24.5 


887 


0.574 


0.398 











are separated by less than 1" were merged, and their V — I 
are computed. 

The completeness magnitudes of the V and / cat- 
alogues are defined to be one half magnitude brighter 
than the peak of the distributions. This corresponds to 
V comp iete — 23.75 and I C om P iete — 22.75. Because some 
chips go deeper than others and because the correlation 
analysis is sensitive to chip-to-chip number density varia- 
tions, we lower the completeness limit to the least sensi- 
tive chip (0.2 magnitude brighter), to which we add the 
chip-to-chip dispersion of (t0.05 mag measured from Fig. 
4. Hence, the final completeness limits are V comp i e t e — 23.5 

and Icomplete — 22.5. 

The colour completenesses in the V and I bands are 
given relative to one another in Fig. || The red galaxy 
completeness limit is determined primarily by how deep 
the V-band data extend. Because the 1^-band catalogue is 
only complete to V ~ 23.5, faint objects redder than V — 
I > 1 will be missed near the limit of the I catalogue, (see 
Fig. ^ and [lO] in section 4.3). This is an important fact to 
be remembered when making colour-selected correlation 
analyses. 

4. Stars and galaxies 

4-.1- Star/galaxy separation 

SExtractor computes a stellarity index for each detected 
object (in the interval 0-1, with 1 for stars, and for galax- 
ies). The stellarity index is determined from a non-linear 



set of equations (Trained Neural Network) (Bcrtin and 
Arnouts, 1996| ). The good seeing of the images (0.7"-0.8") 
allows a robust classification to V < 22 and I < 21. Ac- 
cording to Bcrtin & Arnouts, the algorithm success rate 
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Fig. 6. Ratio of galaxies that are detected in the two 
bands (completeness) versus V and / magnitudes. At faint 
/ magnitudes, the catalogue is biased against red objects 



18 20 22 24 



Fig. 7. SExtractor star/galaxy classification. Stars have 
value 1 and galaxies have value 0. The classification is 
robust up to V < 22 and I < 21 

at these magnitudes is 95% using data with a similar sam- 
pling and a seeing ~ 1". Figure |?] shows the index of the 
~ 30, 000 objects detected both in V and /. For / < 21 or 
V < 22, all objects with an index < 0.9 are classified as 
galaxies. This criterion classifies as galaxies only the ob- 
jects showing a clear evidence of extendedness. For I > 21 
or V > 22 objects with an index< 0.95 are classified as 
galaxies. Because at these faint magnitudes most objects 
have an index< 0.95 (the great majority of objects with a 
stellarity index higher than 0.95 are spurious detections), 
the threshold of 0.95 does not remove the remaining stars 
from the sample. To correct for the fact that most of the 
stars with V > 22 have been misclassified as galaxies, we 
need to apply a correction for the star dilution (see subsec- 
tion on data- induced errors in section ||). We evaluated the 
stellar contaminati on wit h the Galaxy star-count model of 
Bahcall & Soneira ( 1986 ), which is compared on Figure || 
to the number counts of galaxies. The galaxies outnum- 



ber the stars by nearly an order of magnitude where the 
classification algorithm efficiency is less than 95% (vertical 
dotted line in the lower part of the diagrams). 

4-2. Galaxy counts 

The galaxy counts shown in Fig. |s| are in good agreement 
with other measurements. In the I band, We systemati- 
cally measure 20 — 30% more galaxies than Postman et al. 
( |1998| ) up to I < 22. At J = 21.25, we count 5626 galaxies 
per deg 2 , and Postman et al. find 4057. Given the errors 
in the I zero-point calibration (ctaj ~ 0.05 — 0.2), the pos- 
sible difference in the magnitude scale, and the intrinsic 
cosmic variance, we do not consider this difference to be 
significant. 

We model the galaxy counts of Fig. || following the 
method described by Cole, Treyer, & Silk ( |1992| ). They 
give the equations of the volume element, the comoving 
distance and the luminosity distance of the objects for 
three cosmologies (fio = 1> Einstein de Sitter; flo = 0.2 
Open; Q = 0.2, A = 3(1 - Q )H$ Flat, a factor c/H is 
missing in their equa tion o f the volume element for the flat 
A universe). Yoshii ( 1993 ) also details a similar method. 
The luminosity function (LF) is chosen to be similar to 
the CFRS for which LFs have been measu red separately 
for blue and red galaxies ( Lilly et al., 1995 ). Here we ap- 
proximate the LF by its red c omponent. The S chcchtcr 
parameterization (f>(M) is used ( Schechter, 1976 ), 

4>{M)dM = 0.4 In 10 0* 10 x{1+a) dM x exp(-10 x ), 



X = 0A{M* - M), 



(7) 



where M is the absolute magnitude in the V or / band, 
and 4>*, M* and a are the Schechter parameters (Here, 
Mj = -21.5, 4>* = 0.004, a = -1.0). K corrections are de- 
termined using 13-Gy-old ell iptical galaxy template spec- 
tra from the PEGASE atlas ( Fioc and Rocca-Volmcrange 



1997J), between redshifts < z < 2. As already noted by 
many authors, this simple model does not provide a satis- 
factory fit to the number counts at faint and bright mag- 
nitudes simultaneously for any cosmologies in the V-band. 
A better fit would include a more realistic luminosity func- 
tion accounting for both the red and blue galaxy popula- 
tions, and for either a density or a luminosity evolution 
(see section 6.1). As our purpose here is not to model the 
number counts, we limit ourselves to this partial model. 

Figure || shows that our UH8K number counts deviate 
from the predicted counts in non-evolving Einstein-de Sit- 
ter and open universes at I ^ 21 and at V ^ 22. Postman 



et al. ( 1998| ) observed a clear departure from these two cos- 



mological models in their I counts using the no-evolution 
model of Ferguson & Babul ( |1998| ) (F B). All other surveys 



Cowie et al., 1988 


; Driver et al., 1994|; pardner et al. 


1996; Postman et 


al., 1998|; Woods and Fahlman, 1997; 



behaviour, except for the V number counts of Cowie et al. 
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Fig. 8. Compariso n of ou r UH8K galaxy c ounts in V and / with the res ults of Arnouts et al. ( 199£), Co wie et al. ( 1988 ), 
Drivers & Phillips ( |l994|), Gar dner et al. (|l996|) , Postman et al. Jl998|) , and Woods & Fahlman (|l997|) , Mamon ( pL998|) , 
and McCracken et al. ( 2000a ). Stars counts (solid line) are given according to the Bahcall model of the Galaxy. The 
vertical dotted line gives the limiting magnitude below which the star/galaxy classification is reliable (Fig. 0). Galaxy 
number counts are given for three cosmologies in V and / bands (cf text). The counts are in marginal agreement with 
a no-evolution flat A universe. 



( |1988| ), which agree well with the Einstein-de Sitter model 
at V < 22 and deviate only at fainter magnitudes. 

In contrast, Figure |] shows that the non-evolving flat 
A universe model provides a marginal agreement with our 
UH8K galaxy number counts in the / band at our mag- 
nitude limit of / ~ 22.5. In fact, most surveys displayed 
in the right panel of Fig. || agree with this cosmological 
model at I ^ 22, as expected from the results of the CFRS 
( Lilly et al., 1996 ), where giant red galaxies evolve little in 
the redshift interval — 1. In Fig. ||, the deviations from 
the non-evolving flat A universe occur in dataset which 
probe the faintest magnitudes, near / ~ 24. Our sample 
is not deep enough to show a departure from this cosmo- 
logical model. Lidman et al. also find that their / num- 
ber counts (not shown in Fig. p|) are compatible with a 



no-evolution flat A universe (Lidman and Peterson, 1996) 
out to / ~ 21, whereas they deviate from the evolution 
model of Yoshii ( 1993| ) (this model uses a dwarf galaxy 
blue LF, which yields very different K corrections at faint 
magnitudes compared to a no evolution model) . From the 
mentioned surveys, there are indications that evolution 
should be used in the modeling of the / counts at / > 22, 
where significant departures from a non-evolving distribu- 
tion occur. 



There is also marginal agreement of most V number 
counts displayed in the left panel of Fig. [|, including our 
UH8K data, with a non-evolving flat A universe model. 
The deviations from the model occur in a wider range 
of V magnitude, namely at V ^ 22 — 24, depending on 
the data set, and the deviation is greater than for the / 
counts. This may be partly explained by the fact that the 

V number counts are more sensitive to evolution of the 
blue galaxy population than the / counts. 

4-3. Galaxy colours 

Figure ^| shows histograms of galaxies and stars versus 

V — I colours. For galaxies, histograms are given for three 
limiting magnitudes, / < 20, / < 21, and / < 22.5. The 
corresponding median V — I colours are 1.45, 1.57, and 
1.38. Figure |l^ shows V — I colours vs / magnitudes for 
the sample of ~ 30, 000 detected galaxies. A vertical line 
shows the /-band completion limit / = 22.5 and an oblique 
line shows V — I colour limits accessible due the F-band 
completion limit of V = 23.5. There is no visible trend to- 
wards a strong colour evolution of the galaxies with their 
magnitudes, but a robust conclusion is not possible be- 
cause the I sample is depleted in red objects. A natural 
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Table 4. Differential V and /-band galaxy counts (see 
Fig-!)- 



V 


TV 


a at 


/ 


N 


a at 


15.75 


1.4 


1.4 


15.75 


29.8 


6.5 


16.25 


4.3 


2.5 


16.25 


21.3 


5.5 


16.75 


10.2 


3.8 


16.75 


45.4 


8.0 


17.25 


20.3 


5.4 


17.25 


103.6 


12.1 


17.75 


36.2 


7.2 


17.75 


183.1 


16.1 


18.25 


69.6 


10.0 


18.25 


330.8 


21.7 


18.75 


146.4 


14.6 


18.75 


533.9 


27.5 


19.25 


236.3 


18.5 


19.25 


968.4 


37.1 


19.75 


406.0 


24.3 


19.75 


1434.2 


45.1 


20.25 


656.9 


30.9 


20.25 


2355.7 


57.8 


20.75 


964.2 


37.4 


20.75 


3794.2 


73.4 


21.25 


1571.8 


47.7 


21.25 


5626.0 


89.3 


21.75 


2531.7 


60.6 


21.75 


7626.8 


104.0 


22.25 


4058.5 


76.7 


22.25 


10312.0 


121.0 


22.75 


6219.1 


95.0 


22.75 


13972.8 


140.8 


23.25 


9991.9 


120.4 


23.25 


15936.6 


150.4 


23.75 


15665.8 


150.7 


23.75 


12284.4 


132.0 


24.25 


16725.8 


155.7 
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Fig. 10. Galaxy V— I colours vs I magnitude. The vertical 
line indicates the completeness limit / = 22.5. The oblique 
line shows the V — I colours accessible at the completeness 
limit V = 23.5. 



5.1. Estimation of u(9) 




Fig. 9. Histogram of galaxy counts (solid line) from top 
to bottom / < 22.5, / < 21, / < 20 and star counts (dot- 
dashed line) vsV — I colour. The star counts are scaled up 
by an order of magnitude. Only stars brighter than / < 19 
are included. 

although arbitrary value to divide red galaxies from blue 
galaxies is the median of the histogram with / < 22.5 in 
Fig. ||, located near V — I = 1.4. In the rest of the paper 
red galaxies will be those having V — I > 1.4 and blue 
galaxies will be those having V — I < 1.4. 



5. Analysis 



The 2-point angular correlation function lu(8) is calculated 
by generating samples of random points covering the same 
area and having the same number as the galaxy sample. 
We use the estimator W(9) defined by Landy & Szalay 
( |1993| ) (hereafter LS), which has the advantage of reduced 
edge effects and smallest possible variance: 



W(6) = W 



DD - 2DR + RR 
RR 



(8) 



where DD is the number of galaxy-galaxy pairs, D R the 
number of galaxy-random pairs, and RR is the number 
of random-random pair s, all of a g iven angular separa- 
tion 6. Following Roche ( 1996 ; 1999 ), we set a logarithmic 
binning for the separation defined as Alog(0) = 0.2. The 
numerical approach of LS is used to calculate DD, DR 
and RR. If one defines the variables d and x as 



d = 



DD 



G p (6)n(n-l)/2'- 
DR 



G p (0)n 2 ' 
G p (9) =< RR > /[n(n- l)/2], 
then eq. (||) can be re-written as 



W = d - 2 x 



1 



(9) 
(10) 
(11) 

(12) 



G p {9) is the probability of finding two randomly placed 
galaxies separated by a distance 6±d9/2, n is the number 
of real galaxies in the sample, and < RR > is the average 
over N realizations of a random sample of k points, with 
Nk = 100, 000 (k varies between 10 to 100 depending on 
the magnitude interval in which the angular correlation 
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function is measured). LS also define Gt(9), the probabil- 
ity of finding two neighbors both at a distance 9±dd/2 of 
one given object. 



G t {6) =< n t (9) > /[n(n - l)(n - 2)/2] 



(13) 



where < n t (9) > is the average number of unique triplets. 
Gt( 0) is necessary to evaluate the random errors (cf section 
5J3). 



5.2. Modeling of u){6), £(r) and ro 

The canonical parameterizations of the two-point spatial 



correlation function £(r) ( Phillipps et al., 1978) a nd of the 
angular correlation function lo{9) ( Peebles] 1980 ) are 



£(r,z) = 



( r / r ooY 
{1 + z) 3 



and 



lo(0) = A u 0- 



(14) 



where r is the comoving distance, r 00 the correlation 
length at z = 0, 9 the angular separation in radian, is 
the amplitude of angular correlation function, 7 and 6 are 
the slopes (5 = 7—1) and e is a parameter characterizing 
the evolution of clustering with the redshift z. If e > the 
clustering evolves in proper coordinates, if e = the clus- 
tering is constant in proper coordinates hence increases in 
an expanding universe, if e = —1.2 the clustering is con- 
stant in comoving coordinates. The comoving correlation 
length at a redshift z is related to roo by 



r (z) = r 00 (l + z)-( 3+£ -^. 



(15) 



Given eq. (|l4|) and the galaxy redshift distribution 
N(z), one can relate ui(9) and £(r): 



= C r (z peak )i O 1 ^ B(e) 
B(e) - 



c 

#0 



7-1 



ff(*)(l + *)3+« 
2 

N(z)dz' 



dz 



N{z) 

g{z) 
g(z) 

c 



(tt/180)" 



J2MM,z) 



dV 



(1 + Z) 2 y/1 + Q z\ 

(1+z)" 1 



y/Q (i + z y - n 

rK7_i]/2] 
^ r(7/2) " 



+ 1 
3.68 



du! 
(A = 0), 

(A + n = 1), 



(16) 

(17) 
(18) 

(19) 
(20) 

(21) 



Here r<j(z) is the angular diameter distance, dV/duj dz is 
the comoving volume element, both given for three cos- 
mologies in the appendix of Cole, Treyer and Silk ( 1992j ). 
5^. 4>i{M, z) is the luminosity function, whose definition 
might be dependent on different spectral type i evolving 
with z, and Y is the gamma functi on. The value of C is 
given for a typical value of 7 = 1.8 ( Peebles, 1980 ). Equa- 
tions (|16f) to (Ell) allow to make a direct derivation of 



ro(zp e ak) from A u! (Am), where z pea k if the peak of the 
redshift distribution of the galaxies in the sample defined 
by the interval of apparent magnitude Am: 



Au (zp ea k ) 



1/7 



ro(z peak ) = [^fg^j ■ (22) 
(A u (Am) is rewritten as A ul (z pea k))- 

5.3. Estimation of errors 

In the past few years, considerable theoretical efforts have 
been devoted to the calculation of errors in the estimation 
of oj{9). The errors can be divided into two categories; 
(1) the random errors; and (2) the systematic errors due 
to the various observational biases in the data. The sys- 
tematic errors are caused by false detections, star/galaxy 
misidentifications, photometric variations and astrometric 
errors. The random errors are induced by the finite area 
of the survey and depend on the geometry and the size of 
the sample. 

5.3.1. Random errors 

Until recently, a Gaussian approximation for the distribu- 
tion of the galaxies was assumed to calculate the errors on 
w(9). However, the distribution of the galaxies is known to 
depart significantly from a Gaussian distribution on small 
scales. A finer approach would be to include the possible 
correlations of the data into the error analysis. This has 



been done by Bernstein (1994) on the LS estimator W{9). 
Bernstein derives a formal solution to the random errors 
for the case of a hierarchical clustering universe in the lim- 
its n 3> 1 (number of galaxies in the sample), W <C 1, and 
9 <C angular size of the sample: 



AW\ 2 



4 (1 - 2 q 3 + q 4 ) We 



4 
n 
2 



W r (l + 2 q 3 W) 



W 2 



(G- 1 1 



53-1 



1 + W 

W 2 



W' 1 - 1 



(23) 



where the parameters and q^ are measured by 



Gaztanaga (1994). They are related to the hierarchical 
amplitudes S3 and S4 (s n =< oj n > /[< u>2 >]" , 
where w n are the n-point angular correlation functions), 
by <;3 ~ S3/3 and <?4 ~ S4/I6 (Gaztanaga, 1994). S3 and 
S4 have been derived by Gaztanaga (1994f) and Roche & 



Eales (1999) from the APM catalogue ( |Maddox et al. 
1990a| ): s 3 ~ 4 and s 4 ~ 50 at 9 ~ 0.1°, the scale at 
which the angular-correlation function is well measured 
in the APM catalogue ( pVIaddox et al., 1990b| ); therefore 
(73 ~ 1.3 and 94 ~ 3. W r (r for ring) is the 3-point angular 
correlation for triplets of galaxies defined by 2 galaxies at 
a distance in the interval [9,9 + 69] from the 3rd galaxy. In 
principle, W r is marginally greater than W(9) as defined 
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in eq. (B) but W r ~ W is a fair approximation. W ' e ma3 . IS 
the average value of the angular correlation function at the 
largest angular separation of the sample. This quantity is 
not directly accessible because the estimator W is biased 
by finite volume errors at separations similar to the size 
of the survey. To get an estimation of W g max we divide 
our sample into 8 sub-samples (one separation at mid- 
declination and three separations in right ascension) and 
measure the LS estimator W . Then, assuming a power- 
law whose slope is set at the small angles we extrapolate 
the value Wg ma:c for each sub-sample. Finally, Wg max is 
obtained as the average over the 8 sub-samples. Subdiv- 
ing the sample into 8 sub-samples also allows to measure 
the error on W independently at all scales (the largest an- 
gular separation being a quarter of the largest full survey 
separation), by simply taking the variance over the 8 sub- 
samples. The resulting variance largely underestimates the 
variance defined in eq. (p3|). 

5.3.2. Systematic errors 

The cosmic bias or integral constraint: The measurement 
of W$ provides the cosmic bias bw as 



1997) is to apply a star dilution (multiplicative) factor 



D s tar to the parameterized amplitude A u 



b w (e)^(3~4q 3 ~W(9)- 1 ) Wg„ 



(24) 



(Colombi, priv. comm. 1999). This bias can be corrected 
in the calculation of W(9) in the form of an additive fac- 
tor IC. At small angular separations eq. (^ ) simplifies as 
bw — —We maa ./W. This negative bias is very small but 
becomes comparable to W when 6 approaches the angu- 
lar size of the sample. The usual way of correcting for 
the cosmic bias is to fix a slope for the angular correla- 
tion function and to find the constant value of IC which 
minimizes the \ 2 fitting to the data. One may point out 
that the cosmic bias becomes non-negligible only when the 
errors AW/W of the LS estimator given by eq. (|2^) are 
large. The smallest scale at which bw is 10% is 8 > 3', 
which corresponds to the last 2 points in all curves plot- 
ted in Figures 12 and 13. In the interval 17 < I < 21 for 
example, bw — —0.09. However, the random errors at this 
scale are also significant, AW/W ~ 12%, which gives lit- 
tle weight to these points in the least-square fits of W{&). 
We therefore consider that the cosmic bias has negligible 
impact on our reported slope and correlation amplitude, 
and we chose to ignore the cosmic bias in order to avoid 
introducing a prior information on the slope of W(9). 



Misidentified stars: Because stars are uncorrclatcd on the 
sky as shown on Fig. |ll|, the stars fainter than V > 22 and 
/ > 21 which are not removed from the catalogue dilute 
the clustering present in the galaxy correlation, i.e. de- 
crease the am plitude of m(9). The usual way of correcting 
for this bias ( |Postman ct al., 199S : Woods and Fahlman, 



Nobj 



■*Y obj -*Y star 



(25) 



where N O0 j is the number of objects in the galaxy samples, 
and N sta r is the numbe r of stars pred icted by the Bahcall 
model of the Galaxy flBahcall, 1986|) . For V < 22 and 
I < 21, N s tar is given by the classification efficiency of 
95% of SExtractor, namely the upper value is 5% of the 
number of stars detected, so D star < 1.1. The values of 
Dstar for fainter limiting magnitudes are given in Table § 

Masking: Although the SExtractor programme is very 
good at avoiding false detections, it is sometimes tricked 
by the diffraction patterns and large wings of bright stars 
as would be any code using an isophotal threshold for 
detection. Because such false detections are strongly clus- 
tered, they increase the correlation amplitude at scales 
corresponding to the angular size of the false structure. 
To prevent the systematic patterns which could be in- 
troduced by false detections, we define two masks, one 
covering the bad pixels, columns and regions, and one 
covering all stars brighter the V < 15 (from the HST 
Guide Star Catalog). The second mask is made of four 
empirical components, (1) a central disk whose radius is 
defined by an exponential law as a function of magni- 
tude (radius — 4969 exp[— 0.378M ag] pixels), (2) a ver- 
tical rectangle covering the bleeding streak, whose length 
is an exponential function of magnitude (bleed = 1.28 x 
10 7 exp[— 0.8Mag] pixels), (3) and (4) inclined rectangles 
(at 38.5° and 51° of the vertical axis) covering the diffrac- 
tion spikes, also following exponential law of the magni- 
tude (spike = 5292 exp[— 0.315Magl pixels). The masked 
regions are partially visible in Fig. |l| because only bright 
stars show large wings and the low density of objects does 
not permit one to distinguish "true" empty regions from 
masked regions. We apply the same two masks to the ran- 
dom realizations for the evaluation of LS estimator W(0) 
in section 5.1. 



Photometric errors: Two kinds of photometric errors may 
be present: the random errors or photon noise called <jm, 
and the residual calibration errors in the coefficients Ay, 
Aj, ky, ki and the zero-points Vb and Iq. Note that the 
errors given in Table [}] called ctav and cta/ include both 
the random and calibration errors. It is difficult but nec- 
essary to evaluate quantitatively the two kinds of errors 
because they have opposite qualitative effects on the an- 
gular correlation function. 

The random errors can be evaluated from Table |^. It 
is clear that at bright magnitudes, the calibration errors 
dominate while random errors dominate at faint magni- 
tudes. A random error in a galaxy apparent magnitude is 
equivalent to an increase of the possible volume in which 
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that galaxy lies, i.e. it is equivalent to a convolution of 
the de-projected distance interval. Hence, the random er- 
ror erases the clustering present in the sample. It is diffi- 
cult to evaluate directly the decrease in the amplitude of 
uj(8) due to the random errors because it is impossible to 
disentangle it from a real variation of the spatial galaxy 
clustering. Nevertheless, one can follow a simple argument 
to estimate how random errors affect the measurement of 
uj(8). First, assume that the galaxy number counts follow 
a power-law logN = a mag + cnst so the relative error is 
AN /N — (a/ log e)Amag. An extreme case is to consider 
that all the galaxies with a magnitude error superior to 
the magnitude bin for which u>(8) is evaluated (~ 0.5), 
are uncorrelated to the sample actually falling in the bin. 
This is an extreme case because many of the galaxies with 
the large error in magnitude do belong to the bin. In that 
case, these galaxies would have a very similar effect on the 
amplitude of ui(8) as if they were stars. So the multiplica- 
tive diluting factor D ga i of the random photometric error 
would be. 



Dgai = (1 - (a/ log e)Amag)~ 



(26) 



where Amag is the magnitude error and a is the slope of 
the galaxy number count power-law. Taking the best fit 
slope a ~ 0.4 (see Fig. |8| and Table ||), a typical magni- 
tude error of 0.15 (see Table ||) leads to a dilution factor 
of Dg a i — 1-3, of order of the star dilution factor D sta r 
given in Table §. This estimate of D gal is an indicative 
upper limit, and cannot be used to correct for the dilu- 
tion due to random errors in the magnitude of the galaxies 
because the prior condition is that these galaxies are un- 
correlated. This assumption might not be true, and the 
correction would then artificially increase the amplitude 
of the correlation function. 

We now evaluate the calibration error budget. Because 
the largest airmass difference in V is 0.267, and 0.028 in 
/, even a 50% error on Ay and Aj would not produce 
more than AV ~ 0.015 and AI ~7x 10~ 4 . If we assume 
that the chip-to-chip error on fcy is given by the differ- 
ence between the values of the synthetic sequence and the 
measurement of Cuillandre in Table || (this is certainly 
not the case for fc/, because Cuillandre's value is too low) 
then in the most extreme colours (V — I) ~ 4 (only a 
very small fraction of our sample), the resulting magni- 
tude difference would be AV ~ 0.008. The case of fcj is 
not clear, although it seems reasonable to assume that 
the error cannot be more than ten times the error AV, 
so AI < 0.08. The residual systematic errors in the zero- 
points have been evaluated in section 3.2 to be a ~ 0.05. 



Combining all these mentioned sources of calibration er- 
rors, one finds a value of 0.053 in V and an upper value 
of 0.094 in I. 

The residual photometric errors have an opposite ef- 
fect on W(9). Namely, they introduces CCD-to-CCD vari- 
ations in the galaxy number counts wrongly interpreted by 
the correlation analysis as intrinsic clustering on a scale 
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Fig. 11. Angular correlation function of stars brighter 
than I < 21. The error bars show 1 a Poisson noise. The 
result is compatible with a random distribution. 



given by the angular size of the individual CCD chips, 
thus artificially increasing the amplitu de of oj{8) on these 
scales. Geller, de Lapparent & Kurtz ( 1984 ) showed that 
plate-to-plate systematic variations of more than 0.05 mag 
introduced in the Shane- Wirtanen catalogue would pro- 
duce a flattening of the correlation function and an ar- 
tificial break at a scale corresponding to the plate size. 
For that reason, we limit our measurement of the angu- 
lar correlation to 8 < 400" ~ 0.1°, corresponding to the 
smallest dimension of the individual CCD's. We point out 
that the combined effect of random and zero-point errors 
would be to flatten artificially the slope of oj(8), as Geller 
et al. demonstrated. 



Astrometric errors: Small scale astrometric errors are 
likely to become significant when the bin size of u){0) is of 
the order of the error. The trivial method to avoid such 
contamination is to limit the analysis to bins greater than 
the errors, i.e. greater than 0.5" in our case. When one 
builds a mosaic survey combining overlapping patches of 
the sky, systematic astrometric errors may come into play: 
the number of objects may be artificially higher or smaller 
in overlapping regions just because of misidentifications. 
This would induce a similar effect on the amplitude of 
correlation to the zero-point errors at an angular scales of 
the order of the field size. However, a comparative analy- 
sis of the number counts of overlapping regions with other 
parts of the field does not show any significant bias toward 
a higher number of objects in the overlapping regions. We 
therefore consider that our astrometric errors have a negli- 
gible effect on u>(0), and we limit its calculation to 9 > 1". 
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Table 5. Star dilution correction factor D sta r in V and I 
bands 



V 


Dstar 


/ 


-Dstar 


22.25 


1.08 


21.25 


1.23 


22.75 


1.11 


21.75 


1.20 


23.25 


1.10 


22.25 


1.17 


23.75 


1.09 


22.75 


1.15 


24.25 


1.08 


23.25 


1.13 



6. Results 

First, we tested the reliability of our code for measuring 
the LS estimator of oj(9) on the Zwick y cata logue. The re- 
sult is consistent with that of Peebles ( 1974 ) . We also used 
a deep catalogue (courtesy of Roukema, priv. comm. 1999) 
and found good agreement with its uj(0) measurement. 

Figure |ll| shows the result of our angular correlation 
code on the sample of stars brighter than I < 21. We limit 
the calculation of W(9) at 0.1° because of the reasons in- 
voked in section 5. The function is always consistent with 
a random distribution. The weak positive signal might be 
a sign of misclassification of galaxies or simply small over- 
densities. The LS estimator (eq. |l^]) is then measured, 
using the entire sample of galaxies, for six cumulative mag- 
nitude intervals: 17 < / < 20, 17 < / < 20.5, 17 < / < 21, 
17 < / < 21.5, 17 < / < 22, and 17 < / < 22.5, 
and four incremental magnitudes intervals: 20 < / < 21, 
21< / < 21.5, 21.5 < I <22, and 22 < I < 22.5. 

Random errors are computed using eq. (p3|). Table [3] 
gives the values of W(9) and the 3-cr random errors for 
each cumulative magnitude interval, while Table |7] gives 
the best-fitted values of the corresponding W(9) for the 
parameterization of eq. (Mj)- Three values of the ampli- 
tude are given in Table [7|7An/ is the best-fit amplitude 
when both the slope and the amplitude arc fitted. Ayy is 
the best-fitted amplitude when the slope is constrained to 
S = 0.8, and Af*/"" is A®^ corrected for the star dilution 
l^ 8 x D st ar- The amplitudes Aw, 
correspond to the values of the correlation 
function at 1°, the usual chosen reference scale. Table [z] 
shows the results for both the integrated intervals and in- 
cremental intervals of magnitudes. All quoted errors are 
3-cr. For the same I me dian, the values of A^ are consis- 
tent between the integrated intervals and the incremental 
intervals. Because the error bars are larger for the incre- 
mental intervals, we based our analysis on the integrated 
intervals only. Figure [l^ plots the LS W (9) in logarith- 
mic scale for the /-band incremental magnitude intervals, 
along with the best-fit power-laws for an unconstrained 
slope (Table 0). Figure |l3| plots the difference of all W{9) 
with the constrained S = 0.8 power-laws. 

A final check on the code was done by Colombi (priv. 
comm. 1999) using a count-in-cells routine. The results 
showed that the quality of the dataset allows to measure 



factor (I > 21) A% ar 
A™, and Aft" 



higher orders of the distribution of galaxies (skewness & 
kurtosis). The measurements of the angular correlation in 
the interval 17 < / < 21 gave consistent results with the 
LS estimator. 

Here, we only show the correlation function in the / 
band. We also measured W(9) on the F-band map in the 
magnitude intervals 17 < V < 21.5, 17 < V < 22, 17 < 
V < 22.5, 17 < V < 23, and 17 < V < 23.5. The values 
of the amplitude Aw obtained in the V band are very 
similar to the values obtained in I, but the average slope 
S is ~ 0.5 — 0.6 in V, as compared to ~ 0.8 — 0.9 in 
/. Neuschaefer & Windhorst (1995) also measured slopes 
5 ~ 0.5 in the g and r bands and S ~ 0.7 — 0.8 in the i 
band. 

6.1. Variation of slope versus magnitude 

Recent ACDM models predict a decrease of the spatial 
correlation slope 7 for scales < 10/i -1 Mpc from 7 = 1.8 
at z = to 7 ~ 1.6 for z > 1 ( Kauffmann et al., 1999| ). 
This implies a decrease of the slope <5 = 7 — 1 of W(9) at 
small angular scales and faint magnitudes. Observational 
evidences are poorly conclusive. Brainerd et al. ( |1999| ) re- 
port a steepening of the slope on small scales while Cam- 
pos et al., 1995, Neuschaefer and Windhorst, 1995, Infante 
and Pritchet, 1995, and Postman et al. 1998 find the oppo- 
site effect. Other authors find no significant variations to 
the limits of their samples (pouch et al., 1993|; Roche and 
Eales, 1999; Hudon and Lilly, 1996; Woods and Fahlman, 
1997| ) . Figure |13| and Table |7| show no significant flattening 
of slope at faint magnitudes. The slope is compatible with 
5 = 0.8 for all magnitude intervals. Our result is neverthe- 
less consistent with the results of Postman et al. ( 1998 ) 
who find signs of a decrease only in their faintest bins, at 
I > 22, near and beyond our / limit. 

6.2. Variation of Aw with magnitude 

The choice of a galaxy luminosity function to model the 
decrease of Aw with the limiting apparent magnitude is 
crucial because the behavior of Aw is sensitive to both 
the parameterized characteristic absolute magnitude M* 
and the slope a (eq. [Q). We choose to use the luminosity 
function observed in the CFRS (Lilly et al. 1996), derived 
from 591 galaxies in the range < z < 1, keeping in 
mind that such a small sample of galaxies can only provide 
an indication of a general trend. Because both the CFRS 
sample and our sample have been selected in the / band, 
we thus limit the possible biases due to the photometric 
sample selection. 

Lilly et al. divide the CFRS sample into a red popula- 
tion (redder than an Sbc having rest-frame [U — V]ab = 
1.38; [U - V]ab ^[V- I]ab at z ~ 0.5, V AB = V, and 
Iab = I + 0.48) and a blue population (bluer than an 
Sbc). The red galaxies show no density or luminosity evo- 
lution in the range < z < 1, whereas the blue galaxies 
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Table 6. Measures of W{6) and 3-a errors (from eq. [||) in the I band for 17 < 7 < 20 (W 20 ), 17 < I < 20.5 (W20.5), 
17 < I< 21 (W21), 17 < 7 < 21.5 (W21.5), 17 < K 22 (W 22 ), and 17 < I< 22.5 (W 2 2.5)- 



6» 


W20 


3 (T20 


W20.5 


3 (720.5 


W21 


3 (T21 


W21 5 


3 021. 5 


W22 


3 (T22 


W22.5 


3 022.5 


arcsec 












(xlO 


" 3 ) 












1.0 


-753 


538 


-750 


359 


-1068 


242 


-545 


210 


-469 


150 


-377 


110 


1.6 


538 


1041 


-134 


429 


136 


299 


-84.0 


176 


-61.1 


119 


-222 


83.3 


2.5 


1467 


889 


1346 


574 


760 


309 


554 


194 


422 


124 


291 


78.6 


4.0 


430 


389 


650 


301 


570 


204 


380 


123 


334 


86.8 


311 


66.9 


6.3 


767 


356 


561 


216 


446 


143 


392 


105 


287 


66.6 


226 


46.3 


10.0 


464 


212 


359 


134 


291 


90.8 


236 


63.3 


180 


41.6 


146 


29.8 


15.8 


372 


152 


272 


95.0 


267 


75.5 


182 


46.6 


136 


30.3 


121 


23.5 


25.1 


331 


123 


215 


69.7 


171 


48.3 


122 


31.1 


95.3 


20.9 


86.5 


16.6 


39.8 


209 


77.9 


160 


50.3 


127 


35.1 


94.2 


23.5 


71.3 


15.4 


63.4 


12.2 


63.1 


145 


54.1 


112 


35.3 


71.6 


20.7 


56.4 


14.5 


45.3 


10.1 


41.2 


8.06 


100.0 


83.5 


32.4 


70.1 


22.7 


49.7 


14.6 


35.3 


9.41 


30.2 


6.89 


27.5 


5.54 


158.5 


55.8 


22.2 


43.8 


14.9 


26.4 


8.45 


21.5 


6.10 


17.8 


4.36 


17.2 


3.65 


251.2 


30.5 


13.4 


22.0 


8.50 


13.3 


4.95 


10.3 


3.43 


10.7 


2.86 


7.09 


18.2 


398.1 


16.4 


8.38 


3.19 


2.99 


2.30 


1.93 


-1.23 


0.51 


1.61 


0.95 


1.09 


6.37 


631.0 


-2.29 


-1.40 


-4.12 


1.15 


-2.27 


0.68 


-1.80 


0.56 


-0.48 


0.11 


-1.44 


0.40 


1000.0 


-7.00 


-2.18 


-5.14 


1.62 


-2.47 


0.96 


-0.50 


0.16 


-0.64 


0.31 


1.17 


0.56 


1584.9 


-3.04 


1.88 


-1.85 


1.84 


-4.98 


1.26 


-5.00 


0.93 


-3.68 


0.73 


-2.50 


0.53 


2511.9 


-32 


7.91 


-1.58 


2.35 


-9.85 


0.86 


2.27 


1.21 


1.15 


0.33 


-2.56 


0.54 



Table 7. Best-fit values of the amplitude Aw and the slope 5 of W{9) — Aw@~ S ', Ayy is obtained for a fixed slope 
5 = 0.8 and A w a,r is the corresponding star dilution corrected amplitude (see text) in the interval 10" < 8 < 500". All 
quoted errors are 3 a. The first 6 rows of the table show the integrated intervals of magnitudes, the last 4 rows show 
the incremental intervals of magnitudes. 



I Magnitude 


Ngal 


A w 


8 


A w 


A star 
W 


Interval 


Median 




xl0~ 4 




xlO" 4 


xlO" 4 


17-20.0 


19.356 


2166 


40.4 ±61.3 


0.839 ± 0.329 


47.9 ± 19.6 


47.9 


17-20.5 


19.807 


3552 


31.7 ±40.8 


0.830 ± 0.277 


36.1 ± 12.5 


36.1 


17-21.0 


20.281 


5671 


15.6 ±17.8 


0.917 ±0.240 


25.8 ± 8.02 


25.8 


17-21.5 


20.724 


8739 


11.42 ± 11.87 


0.930 ±0.218 


19.95 ± 5.48 


19.95 


17-22.0 


21.166 


13215 


12.6 ± 11.4 


0.861 ±0.190 


16.6 ±3.85 


19.4 


17-22.5 


21.617 


19501 


10.06 ± 7.80 


0.872 ±0.161 


13.8 ± 2.88 


15.9 


20-21.0 


20.610 


3504 


9.35 ± 14.4 


0.981 ± 0.320 


20.0 ± 8.12 


20.0 


21-21.5 


21.270 


3068 


6.67 ± 14.9 


0.947 ± 0.477 


12.36 ±7.05 


15.2 


21.5-22 


21.772 


4476 


10.3 ± 17.6 


0.880 ± 0.372 


14.5 ± 5.82 


17.5 


22-22.5 


22.267 


6283 


11.2 ± 14.0 


0.810 ±0.277 


11.7 ±4.01 


13.7 



show a luminosity evolution of about 1 magnitude in the 
same redshift ra nge. Recent measu rements made from the 
CNOC2 survey ( |Lin et al., 1999| ), on -2000 galaxies in 
the range 0.1 < z < 0.7, confirms the general observa- 
tions of the CFRS, although the proposed interpretation 
is notably different. The CNOC2 analysis separates the lu- 
minosity evolution from the density evolution. Early and 
intermediate (red) galaxies show a small luminosity evo- 
lution in the range 0.1 < z < 0.7 (AM* ~ 0.5), whereas 
late (blue) galaxies show a clear density evolution with al- 
most no luminosity evolution in the same redshift range, 
in apparent contradiction with the CFRS results. 

We choose to adapt the CFRS LF to our sample, and 
we proceed as follows. Galaxies are separated into two 



broad spectro-photometric groups, the E/SO/Sab (called 
red) and the Sbc/Scd/Irr (called blue), using the me- 
dian colour V — I = 1.4 (see section 4.3). The red group 
has a non-evolving luminosity function with parameters 
4>* = 0.0148 h 3 Mpc" 3 , Mj = -21.5 + SloglO/i, and 
a = —0.5 (eq. [Q), and the blue group has a mild-evolving 
luminosity function with parameters </>* = 0.015/i 3 Mpc~ 3 , 

)) + 5 log 10ft, and a = -1.07. 



M; = -21.56 
The factor 1 - 



(1- 

-1z 



-2z 



equals for z = and — 1 for z > 1 
and mimics the observed smooth brightening of Mj with 
redshift. 



We integrate eq. ( |18| ) in the different apparent magni- 
tude intervals listed in Table [| and obtain N(z) for three 
cosmologies. The K corrections are computed from tern- 
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Fig. 12. Plots of Log W(0) spaced by 0.5 dex (symbols; 
see also Table ^) and best-fit curves {Aw and 8 of Ta- 
ble 0) as dotted line. From top to bottom, 17 < I < 20 
(k=5),17 < J < 20.5 (k=4), 17 < I < 21 (k=3), 
17 < 7 < 21.5 (k=2), 17 < 7 < 22 (k=l), 17 < 7 < 22.5 
(k=0) 
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Fig. 13. Plots of the differences of the LS W{6) and the 



best-fit A9, 



Sfl-0.8 



(dotted lines). The curves are spaced by 



1 dcx for sake of clarity. Same magnitude intervals as Fig. 

in 



plates of E (9 Gy) for the red group and Sd (13 Gyr) for 
the bl ue group fro m the PEGASE atlas ( Fioc and Rocca- 
Volmei ange, 1997 ) and are shown in Fig. 14. The choice 
of a given atlas is not benign, as Galaz ( 1998 ) shows that 
K corrections can vary by 50% at z ~ 1 when comparing 
the P EGASE atlas with the GISSEL atlas QBruzual and 
Chariot, 1993), leading to significant differences in N(z). 

Figure 15 shows N(z) for red and blue objects having 
7 < 22.5, for three cosmologies: Einstein-de Sitter fi = 1 



as a solid line, open universe f2 = 0.2 as a dashed line, 
and flat A f2o = 0.2 universe as a dotted line. One can note 
that red and blue galaxies show very different redshift dis- 
tributions as expected from the different LFs. This should 
be kept in mind when we compare the different evolution 
of Aw with magnitude for the red and blue samples. The 
resulting N(z) used to model our sample is the sum of 
the N(z) for the red and blue object distributions respec- 
tively. As our 2 colour samples suffer from differential in- 
completeness (see section 3.3), we normalized the relative 
number of blue and red objects to the observed ratio in the 
CFRS survey at the corresponding limiting magnitude. 

Figures |l6|, [l7|, and 18 show the decrease of the ampli- 
tude of the correlation function corrected for star dilution 
A w ar for the median 7 magnitude of the integrated and 
incremental intervals (Table ^|); a fixed slope of 8 = 0.8 
is used to measure the reference scale, which is taken a 1 
degree. The data-points for the incremental intervals show 
larger error bars, but are consistent with the amplitudes 
for the integrated intervals. The measurements of Postman 
et al. ( 1998] ) and Lidman & Peterson ( 1996| ) are shown for 
comparison. We also plot the expected curves for the three 
universes mentioned above, for each of the three values of 
the clustering parameter e = —1.2 in Fig. 16, e = in Fig. 
|l7| , and e = 0.8 in Fig. |lj). For each value of e, the theoret- 
ical curves are calculated with 8 = —0.8 and the best-fit 
spatial correlation length roo at z = 0, using eq. (15) and 
(22). These values along with the corresponding \ 2 of the 
fit are listed in Table || (only the integrated intervals of 
magnitude are used for the fitting, 6 data-points). Note 
that in Fig. [l(], [It], and [l8], different values of roo shift the 
theoretical curves by constant values in the Y direction. 

Several conclusions can be drawn from Fig. [ll| [l7], |l8[ 
and Table ||. We can always find a set f^o, e and r$ which 
fits our data. A universes with f2 m = 0.2 and A = 0.8 
slightly favour null e, in good agreement with the results 
of Baugh et al. (1999) for semi-analytical models of bi- 
ased galaxy formation. In contrast, Table || shows that 
Einstein-de Sitter universes favour positive e (recall that 
e = —1.2 means no evolution in comoving coordinates, and 
e = no evoluti on in physical coordinates) , as obtained by 
Hudon & Lilly ( 1996 ) (e = 0.8) for th e hierarchical cluster- 
ing CDM model of Davis et al. ( 1985|) in Einstein-de Sitter 
Universes. Third, positive values of evolution parameter e 
give better fits to our observations (Fig. |l8|) than nega- 
tive e. Moreover, comparison of the roo listed in Tabic || 
with the local values of roo — 4 — 8ft." 1 Mpc at z ~ also 
suggest null to mild clustering evolution. 

Finally, Table ^ gives the peak redshift z pea k derived 
from the redshift distributions N(z) (using the CFRS 
model luminosity function) for all magnitude intervals and 
for the three cosmologies. Using eq. the values of 

fo(zpeak) are computed for all z pea k, and can be com- 
pared to other results. At z ~ 0.5, we measure a value 
of ro(zp e ak) in the range 3.4 — 3.7 h^ 1 Mpc, depending 
on the cosmological model and the evolution index e. The 
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Fig. 14. K corrections for early-type and late-type galax- 
ies in V & / band. We use the simplest hypothesis of 
non-evolving galaxy spectra over the range < z < 2. 



typical 1-cr uncertainty on the values of ro(z pea k) listed 
in Table ^| is ~ 0.35 h~ x Mpc. The additional error origi- 
nating from the uncertainty in the cosmology and in e is 
estimated from Tables || and || to be ~ 0.35 hr x Mpc. By 
adding these errors in quadrature, we find an estimated 
total error in the correlation length ro of ~ 0.5 hr 1 Mpc. 
This more realistic error can also be applied to the values 
of roo given in Table |[ 

If we assume a flat A = 0.8 cosmology, Table || 
gives ro(zp e ak) — 3.5 ± 0.5 ft-" 1 Mpc at the peak redshift 
Zpeak ~ 0.58 of the / < 22.5 survey. The corresponding 
value at z pea k ~ 0.50 is r^{z pea k) ~ 3.7 ± 0.5 hr x Mpc. 
Within the error bar, this result is in agreement with most 
other angular correlation measurements at z pea k — 0.5 



(Postman et al., 1998; Hudon and Lilly, 1996; Roche and 



Eales, L999; Woods and Fahlman, 1997) 



If we compare with the direct spatial measurements, 
our result is closer to the CNOC2 results than those from 
the CFRS: Carlberg et al. fi nd that the 2300 b right galax- 
ies in the CNOC2 survey ( |Yee et al., 1996| ) show little 
clustering evolution in the range 0.03 < z < 0.65 with 
rg ~ 3.5 — 4.5 h~ x Mpc at z ~ 0.5 , depending on the 



cosmological parameters (Carlberg et al., 2000); whereas 
Le Fevre et al. measure ro — 1.5 ft -1 Mpc at z ~ 0.5 for 
591 I-selected galaxies, implying a strong evolution from 
the local values (r ~ 4 — 8/i _1 Mpc at z ~ 0) (Lilly et al 
1995[)~|As mentioned above, the small correlation length 
found in the CFRS survey may be due to the cosmic vari- 
ance which affects small area surveys. On the other hand, 
neither the CNOC2 survey nor the CFRS survey detect 
an evolution in the slope 7 ~ 1.8, in good agreement with 
our UH8K data. 



Fig. 15. Redshift distributions N(z) for red and blue ob- 
jects assuming CFRS-like LF. The limiting magnitude is 
I < 22.5. The peaks of the red+blue distributions are 
normalized to 1000. The relative number of blue and red 
objects are normalized to the observed ratio in the CFRS 
survey. 
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Fig. 16. The amplitude of correlation A^ ar for the UH8K 
data is plotted against the median / magnitude for the 
same intervals as in Fig. |lj and [l3] (filled symbols) and is 
compared to the results from other groups. The four filled 
squares with large error bars are the measurements for the 
incremental intervals, and the diamonds show the values 
for the integrated intervals (see also Table [7]) . Theoretical 
curves are calculated for e = —1.2 and the corresponding 
best-fit roo (listed in Table ||) obtained with 5 — —0.8, 
and for three universes: Einstein-de Sitter as solid lines, 
flo = 0.2, A = as dashed lines and £Iq — 0.2, A flat uni- 
verse as dot-dashed lines. We use the bimodal CFRS lumi- 
nosity function with mild luminosity evolution described 
in subsection 6.2. 
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Table 9. Rcdshift 



of the peak of the redshift dis- 
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median 



Fig. 17. Same as Fig. 16, but with e = for the theoretical 
curves. The corresponding best-fit roo are given in Table 
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Fig. 18. Same as Fig. |T^, but with e = 0.8 for the theo- 
retical curves. The corresponding best-fit roo are given in 
Table § 

Table 8. Best-fit values of the spatial correlation length 
roo for fixed e and three cosmologies. 



TOO 


c 


Cosmology 


x 2 


h- 1 Mpc 






A 




3.53 ±0.28 


-1.2 


i 





0.0115 


3.49 ± 0.41 


-1.2 


0.2 





0.0125 
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tribution N(z) and the corresponding best-fit correlation 
lengths ro(zp ea k) for different magnitude intervals and cos- 
mologies. Here e = 0.8 
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3.16 
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6.3. Variation of Aw with galaxy colour 

We also calculate the variations of Aw with depth for the 
blue and red sub-samples of our UH8K data, which are 
defined in subsection 4.3. A map of the 8986 red galaxies 
(V - I > 1.4) with / < 22 is shown in Figure [l9] (top 
panel); the 7259 blue galaxies (V — I < 1.4) to the same 
limiting magnitude arc plotted in the bottom panel of Fig- 
ure |l9|. One can see that red objects are more clustered 
than blue objects, as partly reflected by the morphology 



density relationship (Dressier, 1980). Note that the surface 
density of the blue galaxies increases by a factor of ~ 2 
with increasing right ascension. This is probably caused by 
a systematic drift in the photometric zero-point along the 
survey R.A. direction, which was not completely removed 
by the matching of the magnitudes in the CCD overlaps 
(see subsection 3.2). We then perform an angular correla- 
tion analysis on each colour-selected sample. For the blue 
sample, we introduce in the random simulations, prior to 
masking, the mean R.A. gradient measured from the data. 

Figure |2(] shows the decrease of Aw for the red and 
blue galaxies to a limiting magnitude / < 22.5 (see section 
3.3 & 4.3). In principle, the separation should be done on 
rest-frame colours and not on observed colours, but this 
requires prior knowledge of the redshift of the objects. The 
effect of using observed colours is to decrease the resulting 
angular correlation because galaxies of different intrinsic 
colours at different distances are mixed together. We use 
the same angular binning as in the previous analyses for 
consistency (subsections 6.1 & 6.2). Table [lO] gives the 
median magnitude of the red and blue samples. The last 
interval given in Table [l(] includes all 11,483 red galaxies 
(I med = 21.670) and 20,743 blue galaxies (I med = 22.902) 
fainter than I > 17 (to V ~ 24 and I ~ 23). 

Note that the blue sample reaches fainter I magnitudes 
than the red sample. As seen in the section 3.3, this is a 
selection effect due to the fact that only galaxies detected 
in the 2 bands are shown in the sample, and the blue sam- 
ple is complete to V ~ 24. Thus at / = 22.5, the galaxies 
redder than V — I > 1.5 (in fact most of the red sample) 
show sparse sampling. Each interval in Table [lO] is respec- 
tively complete for galaxies having V — / < 3, 2.5, 2, 1.5. 
We emphasize that the last two intervals of the red sam- 
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Table 10. Median magnitude Imedian of the red and blue 
samples for the different cumulative / magnitude intervals. 
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Fig. 19. Map of 7069 red galaxies with observed V — I > 
1.4 (top) and 4589 blue galaxies with observed V — I < 1.4 
(bottom) to I < 22. 

pie (17 < I < 22 and 17 < I < 22.5) are probably too 
incomplete to be considered as fair samples. The incom- 
pleteness induces two competing effect. On the one hand, 
it increases Aw, because the number of objects is artifi- 
cially small and the resulting correlation is higher; on the 
other hand, it makes the median magnitude brighter (cf 
Table |l0|) , thus inducing a steepening of the slope of the 
Aw decrease. In Fig. the obviously erroneous offset of 
last point for the blue sample, corresponding to the full 
magnitude range I > 17 and a median I magnitude of 
~ 22.9, also illustrates the error on the measurements of 
Aw one can expect in extreme cases where magnitudes 
are poorly defined and objects are under-sampled. 

Keeping in mind the warning of the previous para- 
graphs, and the fact that the error bars prevent us from 
drawing strong conclusions, it is striking to see the differ- 
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Fig. 20. Same as Fig. Eq. The galaxies of our sample are 
divided into red galaxies with V — I > 1.4 (filled octagons) 
and blue galaxies with V — I < 1.4 (filled tri angle s). The 
results for the full sample of Postman et al. ( [1998 ) (open 
squares) are shown for reference. Best-fit models are shown 
for a A flat universe with e = 0.8 for red galaxies using 
CFRS red LF (solid line) and blue galaxies using CFRS 
blue LF with mild evolution (dashed line). 



ent amplitude of Aw for the two colour samples shown in 
Fig. |2(]. The larger Aw for red galaxies than for blue galax- 
ies for all Imedian is in agreement with the well-known 
higher clustering amplitude of early-type galaxies (see for 
recent spatial measurements in the Stromlo-APM survey 
by Loveday et al., 1995), which is partly reflected by the 
morphology density relationship (Dressier, 1980). Using 
the first three data-points of the red sample and blue sam- 
ples in Fig. |2(| we calculate that the two distributions 
are different at a 4-cr level. The corresponding clustering 
amplitudes are ro(z pea k) — 5.3 ± 0.5 hr x Mpc for the red 
sample and ro(z pea k) — 1.9±0.9 h^ 1 Mpc for the blue sam- 
ple; both are based on a A universe with e = 0.8. For the 
purpose of comparing these values of ro(z pea k), the men- 
tioned errors assume that the correlation functions for the 
red and blue galaxies have the same slope S ~ 0.8. This 
may not be true (see Loveday et al. 1995), but our UH8K 
sample does not allow to address this issue. We also ig- 
nore the systematic effect of the cosmology which would 
affect the values of ro(z pea k) for the two populations in the 
same direction. The resulting difference in the measured 
To{z pea k) for the red and blue populations are at the 3-ct 
level. 

Note that if the clustering of the red and blue galax- 
ies were identical, the larger average distance for the red 
galaxies compared to the blue galaxies (see Fig. 15) would 
yield lower Aw for the red galaxies than for the blue galax- 
ies. The opposite effect is observed. Moreover, if the larger 
number of red galaxies (in the first 4 magnitude intervals 
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in Fig. g0|) in the red sample (8986) compared to the blue 
sample (7259, see Table [l0]) is caused by large random er- 
rors at the limit of the catalogues, as seen in section 5.3.2, 
these would tend to dilute the Aw for the red sample. The 
detected increased clustering strength for the red sample 
over the blue sample is therefore a lower limit on the am- 
plitude difference with colour. The apparent flattening of 
Aw with Imedian for red galaxies in Fig. ^ may be due 
to the incompleteness of red objects at faint magnitudes 
as discussed above. 

The difference in clustering amplitudes which we mea- 
sure for our red and bl ue sam ples agrees with observ ations 
by Neuschaefer e t al. (1995 ), Lidman & Peterson ( 1996 ) 
and Roche et al. ( 1996| ). Neuschaefer et al. find that disk- 
dominated galaxies (blue V — I) have marginally lower 
Aw than bulge-dominated galaxies (red V — I) using HST 
multi-colour fields. Similarly, Roche et al. observe a 3<t dif- 
ference in Aw for a sample divided into objects bluer or 
redder than b ~ r — 1.64. Lidman & Peterson see a weak 
difference between two samples separated by V — I = 1.5. 
Other authors don't see any difference between blue and 



red-selected samples, such as Woods and Fahlman (1997) 
for a separation of V — I = 1.3, Brainerd et al. ( 1995] ) for 
(ff — r ) = 0-3, Le Fevre et al. ( |1996| ) in the spatial correla- 
tion length of the CFRS for r est-frame (U-R)ab = 1-38, 
and Infante & Pritchet ( |1995| ) for (J—F) = 1. In all these 
cases excepting Infante & Pritchet, who used photographic 
plates, both the number of galaxies and the angular scale 
of the surveys are small. It might be possible that in these 
surveys cosmic variance hides a weak signal. 



7. Discussion 

7.1. Limber's formula 

Limber's formula (written here as eq. [15] and [16]) re- 
lating lu(9) and £(r) is strongly dependent on the shape 
of the redshift distribution N(z) which depends on the 
characteristic absolute magnitude M* and slope a of the 
luminosity functions of the different types of galaxies (eq. 
|Ls| ) . Locally, these parameters cover a wide range of val- 
ues with regard to the environment and the morphological 



types of the galaxies ( 


Binggeli et al., 1988; 


Lin et al., 1999; 


Marzke et al., 199S; 


Bromley et al., 199£ 


; 


Galaz and de 



Lapparent, 2000). No reasons lead us to think that galax- 
ies would show less diversity at higher redshifts. Recent 
observations (in UV, IR, X, and 7) of space-borne obser- 
vatories do provide evidence in this direction as well as 
the discordant measurements of the LF by the CFRS and 
CNOC2. Hence, a prior knowledge of the detailed luminos- 
ity functions is necessary to go beyond a phenomenological 
description of the clustering of galaxies. 

Our computed ro is credible only if the correct lumi- 
nosity functions have been used in Limber's formula: this 
would guarantee that the modeled N(z) is close to the 
true redshift number distributions for each galaxy type. 



Segregating galaxies into red and blue samples based on 
observed colour, as we do here, is also a crude first step, 
and should rather be performed using intrinsic colour or, 
even better, spectral type (cf section 6.3). These in turn 
would require knowledge of the galaxy redshifts. Approx- 
imate redshifts can also be obtained along with spectral 
type for multi-band photometric surveys using photomet- 
ric redshift techniques ( Koo, 19'99] ). Because none of the 



required functions and distributions are available for our 
UH8K sample, we emphasize that the reported results can 
only be taken as phenomenological, and all the comments 
on the deduced clustering must be taken with caution. 

7.2. The cosmological parameters 

Most authors take for granted that different cosmological 
parameters only lead to minor differences in the evolution 
of clustering, compared to the effects due to the uncertain 
luminosity function. Figure [16], to do show that for a 
given luminosity function different cosmological parame- 
ters induce different values of ro and e. From our models, 
ro differs by more than 15% between the A fiat and Open 
universes. These differences cannot yet be distinguished 
with the present data (up to z ~ 1). The dispersion be- 
tween the different observations (Fig. ^l| ) precludes any 
derivation of the cosmological parameters. 

Given a CFRS luminosity function, the A flat universe 
gives the better fits to galaxy number counts, and cluster- 
ing evolution of our sample (e = — 0.8). This is in agree- 
ment with a current (though controversial) interpretation 



of recent type la supern ovae results (Schmidt et al., 1998; 
Perlmuttcr et al., 1999|) . Notwithstanding the numerous 



modern methods to measure the cosmological parameters, 
the present analysis shows that future surveys containing 
10 6 galaxies with known luminosity functions per galaxy 
type and redshift interval to z ~ 1 will be required to pro- 
vide good constraints on cosmological parameters using 
this technique. 

7.3. The evolution of clustering 

Figure |l] shows the decrease of the amplitude of the an- 
gular correlation of our sample compared to most of the 
other recent measurements made in the / band (Brain- 



3rd and Smail, 1998 




Campos et al., 1995; 


Woods and 


Fahlman, 1997; 


Postman et al., 1998; 


Lidman and Peter- 


son, 1996; 


Neuschaefer and Windhorst, 1995; 


McCracken 



3t al., 2000b ). We applied a magnitude translation be- 
tween Neuschaefer & Windhorst i magnitude and our / 
magnitude of I = i — 0.7. 

In the magnitude range 1=20-22, our results are in 
good agreement with most of these results except those 
of Campos et al, which have the highest amplitude, and 
on the low side, those of Lidman & Peterson and of Mc- 
Cracken et al., 3 times lower in amplitude. Postman et 
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Fig. 21. Same as Fig. [L6|. Our measurements are compared 
with recent measurements found in the literature. 



al., Woods et al., and Neuschaefer & Windhorst obtain 
intermediate values at Imedian > 22. 

Postman et al. observe a flattening of Aw for scales 
> 1' and / > 21. The measurements of Brainerd & Smail 
extend the flattening of Aw observed by Postman et al. to 
I r*j 24. We also find a possible flattening in the decrease 
of Aw beyond I ~ 21. The most probable explanation 
for such high dispersion of about a factor ten between 
the different measurements of Aw is the dependence of 
Aw on the square of the number density. Cosmic variance 
may account for most of the discrepancies. The rest may 
be attributable to systematic errors of the estimators due 
to different spatial or magnitude samplings. 

As pointed out by Neuschaefer & Windhorst ( 1995| ), 
a flattening of the slope S (see eq. jl4j]) with the redshift 
or with apparent magnitude would lower the theoretical 
curve for Aw derived with Limber's formula (eq. |l6[]). 
Hence, smaller Aw would still be compatible with a small 
e, or a high Qq. In other words, clustering would grow 
faster at smaller scales than at larger scales. The flatten- 
ing of the spatial correlation function is predicted by N- 
body simulations (Davis 1985; Kauffmann 1999a; 1999b). 
Neuschaefer & Windhorst parameterize the flattening of 
the slope j(z) with redshift z as: j(z) — 1. 75(1 + z)~ c , 
where C ~ 0.2 ± 0.2. Postman et al. (|l998j ) derive 
C = 0.35 ±0.10 with a slightly different parameterization: 
7 (z) = 1.8(1 + z)~ c . 



8. Conclusions 



o The amplitude of the angular correlation function of 
the complete sample decreases monotonically through the 
entire range of magnitude intervals. 

o The flattening in the decrease of the amplitude, ob- 
served by Postman et al., is marginally confirmed by our 
analysis. 

o The best model to fit the evolution of the amplitude 
of our sample is the combination of the CFRS luminos- 
ity function with mild luminosity evolution of late-type 
galaxies and no evolution of early- type galaxies, a A flat 
universe, a clustering evolution with e > 0, and a comoving 
correlation length of tq ~ 3.7 ± 0.5 hr 1 Mpc at z ~ 0.50. 
This in agreement with the local measurements of tq with 
the clustering evolution predicted by CDM hierarchical 
clustering models. 

o Red-selected galaxies show higher amplitudes of cor- 
relation than blue selected galaxies. 

The deep multi-band photometric surveys which are 
in preparation should determine whether these observa- 
tional results on the evolution of clustering are due to 
an inadequate definition of the luminosity functions of 
the different types of galaxies or whether the actual clus- 
tering differences reflect different formation histories of 
disk-dominated vs bulge-dominated galaxies. Ideally, one 
would like to measure the spatial two-point correlation 
function for each galaxy type, and for different redshift in- 
tervals. The luminosity functions and their evolution with 
redshift must be measured accordingly in order to closely 
model the observed redshift distribution of the sources. 
Application of photometric redshift techniques (Arnouts, 
2000 priv. comm.) to the deep extension of the EIS survey 
( |da Costa fc Rcnzini, 1999| ) (http: // www.eso.org /science 
/eis ) should provide new constraints on these functions. 
Another deep survey which will also allow to address these 
issues is the LZT survey ( Hickson et al., 1998| ), which will 
provide accurate rcdshifts to a z ~ 0.05 and reliable spec- 
tral types for ~ 10 6 galaxies to z ~ 1. This will allow a 
more detailed study of the evolution of Aw ■ Note that the 
measured evolution of the clustering amplitude with red- 
shift in these surveys might also provide useful constraints 
on the cosmological parameters. 
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In this paper, we present measurements of the angular 
correlation function for a sample of ~ 20,000 galaxies to 
I < 22.5, and V < 23.5 observed with the CFHT UH8K 
mosaic CCD camera over a contiguous area of ~ 30' x 90'. 
The main conclusions are the following: 
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